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We describe the fundamental equations for description of electron-photon in- 
teractions in insulating nanostructures. Since the theory for the nanostructures 
are based on the theory for the bulk crystals, we shall describe elements of the 
theory for the bulk crystals |l|-[|| in the first three subsections. It is important in 
discussions on optical phenomena that both the envelope function and the Wannier 
functions (or, equivalently, cell-periodic parts of the Bloch functions) should be 
taken into account. The theory for the nanostructure is discussed in the fourth 
subsection. In the final subsection, we briefly mention quantum optics in nanos- 
tructures. 

4.2.1 Electron and Hole Operators in Insulating Solids 

To describe quantum states in solids, it is customary to start from a mean-field 
approximation to electrons in a perfect crystal. Then, according to the Bloch 
theorem, a single-body state of an electron in a solid crystal of volume V takes 
the following form: 

<Mr) = ^^(r), (4.2.1) 
yv 



where b labels bands (conduction band, heavy-hole band, and so on), and Ubk is 
the cell-periodic function which is periodic from cell to cell. The spin index is 
omitted for the sake of simplicity. The normalization is taken as 

/ d 3 r ^ k (r)^(r) = <W<W, (4.2.2) 
Jv 

- f d 3 p < k (p) Myk (p) - S W , (4-2.3) 

where p denotes the coordinate vector in a unit cell, whose volume is vq, and Jd 3 p 
is over the unit cell. 

We can construct the Wannier function which is localized around cell I as 



w b (r - R e ) = J 77 E e _ik ' R< <Mr), (4.2.4) 




where is the center position of the cell, and J^' denotes the summation over 
the first Brillouin zone. This function is normalized as 

/ d 3 r w* b (r-R e )w b ,(r-R e ) = 8 W 6 U >, (4.2.5) 
Jv 



1 A slightly shortened version has been published as section 4.2 of Mesoscopic Physics and 
Electronics, T. Ando et al. eds., Springer (1998), pp. 140-155. 
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and the inverse transformation is 

4>6k(r) = J^5> ik - R ^ b (r-R,). (4.2.6) 

Note that both the Bloch functions and the Wannier functions form complete sets: 

EE'^Wfe(r') = *(r-r'), (4.2.7) 
b k 

53 53 ^*( r ~ ( r ' - = 5 ( r - r ')- ( 4 - 2 - 8 ) 

b e 

To describe many-body states, we introduce the electron field operator tp(r) which 
satisfies 

{i>(r),ft(r>)} = S(r-r'), {^(r)^(r')} = 0. (4.2.9) 
We can expand ip(r) as 

^( r ) = 12 HV6k(r)c 6k = ^ J] Wb (r - R e )c M . (4.2.10) 

6 k b I 

We can easily show that c&k and satisfy 

{c&k,4'k'} = <W<W, {c&k,C(,'k'} = 0, (4.2.11) 

{c b e,cl>£>} = s bb>5 u >, {c M ,c b >e>} = 0, (4.2.12) 

and we can interpret c£ k and as the creation operators of an electron whose 
wavefunction is (f>bk (r) and wj,(r — Rf), respectively. 

We here consider a non-doped semiconductor or dielectric material. It is cus- 
tomary to define the creation operator of a hole by 

^IkT = Ctj-k.1 when v is a valence band, (4.2.13) 

where only in this equation we recover the spin index for the reader's convenience. 
The ground state is the vacuum of the holes and the conduction electrons: 

c c k\g) = h v k\g) = for all c, v and k, (4.2.14) 

where c and v stand for conduction and valence bands, respectively. We hereafter 
call a conduction electron simply an electron. 

It is convenient to decompose the fundamental field tp into the electron and 
hole parts, ip e and ^7,, respectively, as 

^(r)=&(r)+$(r), (4.2.15) 
V- e (r) = 53 53'</>ck(r)c ck , (4.2.16) 

c k 
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€( r ) = E E'^(r)c„ k = E E'^(r)^_ k . (4.2.17) 



t> k f k 

The anti-commutation relations are evaluated as 



{Mr), ^l(r')} = E E W)<^k(r') = <5e(r - r') = 5* e (r' - r), (4.2.18) 

c k 

{^(r), ^t(r')} = E E'^k(r)^k(r') = 5 h (r - r') = ^(r' - r), (4.2.19) 

v k 

and other anti-commutators are zero. Note that 

<5(r - r') = 6 e (r - r') + 6* h {v - r'). (4.2.20) 

Hence, S c and S v work as projection operators onto subspaces which are spanned 
by </> c k(r)'s and 0„k(r)'s, respectively. It is only in each subspace that they work as 
delta functions and ipl(r) and V>^(r) can be interpreted as operators which create 
an electron and a hole, respectively, at the position r. One should always keep 
this limited meaning of S e , 8h, ip e and iph m niind to avoid inconsistencies. 

4.2.2 Effective-Mass Approximation 

A single-electron state can be obtained by operating cj k on \g). It is an eigenstate 
of the mean-field Hamiltonian H m f of the perfect crystal: 

W m /4lff>=£c(k)4|5>, (4.2.21) 

where we have taken ii m f\g) — 0, and e c (k) is called the band energy. A single- 
hole state, on the other hand, is obtained by operating on \g), and its energy 
is -e„(k): 

n m fht k \g) = -s v (k)hl k \g), (4.2.22) 

Suppose that e c (k) takes the minimum value at k = k c , and e v (k) is maximum 
at k = k v . The minimum energy cost for creating an electron in band c and a 
hole in band v is 

e c (k c ) - e v (k v ) = E Gcv , (4.2.23) 

which is called the energy gap between bands c and v. In particular, the energy gap 
between the bottom (lowest) conduction band c m and the top (highest) valence 
band v m is simply called the energy gap, and is denoted by Eq. For simplicity, 
we hereafter assume that this gap is located at k = 0, i.e., k Cm = k Vm = 0: 

£c m (0)-£„ m (0) =E G . (4.2.24) 

We are interested in states with small k in bands c m and v m . We hereafter write 
c m and v m simply as c and v. 

When an additional single-body potential V\, which may be due to impuri- 
ties, external potential, and so on, are present the additional term appears in the 
Hamiltonian; 

H = H mf + f d 3 r ^(r^r^r) = H m f + H\. (4.2.25) 
Jv 
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In this case, c^ k |g)'s and hl k \g)'s are no longer eigenstates, and we must take their 
linear combinations to obtain true eigenstates. The linear combination must be 
taken over different bands and different k 's, in general. However, the band mixing 
may be negligible if (i) the spatial variation of V\ is slow (i.e., the characteristic 
length of the variation 3> lattice constants) and if (ii) the band of interest is well 
separated around k = from other bands. (This condition is usually satisfied for 
the conduction band. For the valence band, however, the degeneracy of heavy- 
and light-hole bands causes band mixing in many semiconductors. This point will 
be discussed later.) 

In such a case, a single-electron eigenstate, 

H\eq) = e er} \er,), (4.2.26) 

which is labeled by some set of quantum numbers rj, takes the following form: 

\eq) = [ d 3 r Vcn (r)^ (r)\g) = [ d 3 r Vcr] {v)^ e {v)\g) , (4.2.27) 
Jv Jv 

<p cv (r) = ^ F C7) (k)0 ck (r) = V^Y, Fc V (Ke)w c (r - R e ), (4.2.28) 
k i 

where 

F cn (r) ee -^£> r F C7) (k). (4.2.29) 

k 

The probability amplitude of finding an electron at the position r in the subspace 
which is spanned by </> c k(r)'s is evaluated as 

(<# e (r)M) = ^(r). (4.2.30) 

Therefore, we can interpret ip e (r) as the single-body wavefunction of the electron 
in the subspace. It is normalized as 

1 = (ev\ev) = I d 3 r |^ C)) (r)| 2 = y'\F cv (k)\ 2 = v V \F CT ,(R e )\ 2 = / d\ \F cri {v)\ 2 
Jv k t Jv 

(4.2.31) 

In a similar manner, a single-hole eigenstate takes the following form: 

H\h V ) = e hv \h V ) (4.2.32) 

\h V ) = [ d 3 r ri n (r)4,(r)\g) = f d 3 r ^(r)^(r)| 5 ), (4.2.33) 
Jv Jv 

<p vv (r) = ' F vr ,(k)4> vk (r) = ^J2F vr ,(R e )w v (r - R £ ), (4.2.34) 

k i 

F vv (r) ee JL^'e^i^k), (4.2.35) 
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where fj is the "hole-conjugate" of rj. That is, we dehne fj in such a way that it can 
be naturally interpreted as a set of quantum numbers of the hole: for a hydrogen- 
like state, for example, r\ = (n,l,m,cr), where (n,£,m) are the orbital quantum 
numbers and a denotes spin, we then take fj = (n,£,—m,—cr). The probability 
amplitude of finding a hole at the position r in the subspace which is spanned by 
0„k(r)'s is given by 

(g\Mr)\hv) =^(r). (4.2.36) 

Therefore, we can interpret (^*j=(r) as the single-body wavefunction of the hole 
which has the set of quantum numbers fj. It is normalized as 



1 = {hr}\hrf) = / d 3 r \<p vv (r)\ 



(4.2.37) 

Under our assumptions (i) and (ii), Fb n (r) (b — c,v) becomes a slowly varying 
function of r, and hence called the "envelope function." We note that if a function 
A(r) is slowly varying and if a function B(r) is a cell-periodic function, then 

/ d 3 r A(r)B(r) ~ V A(R e ) [ d 3 p B{p) ~ / d 3 rA(r) ■ — [ d 3 p B(p) 

JV I Jv Uv J l v Jv 

(4.2.38) 

where J d 3 p means the integral over a unit cell. Using this relation and the 
assumptions (i) and (ii), we can derive the approximate eigenvalue equations for 
the envelope functions as 



[e c (-iV) + V x {y)) F cv {r) = e er) F cr ,(r), 

[-e v (-iV) - Ui(r)] F ,(r) = e hn F vv {r). 
In particular, when the bands have parabolic dispersions; 

e c (k) = e c (0) + 
e v (k)=e v {0)- 



2m e 

h 2 k 2 

2mh 



then we obtain the "Schrodinger equations" for the envelope functions: 



2m. 



-V 2 +£ C (0) + V x (r) 



Fc?y(r) — ^cri^cri (r) , 



2m, 



-V 2 -e^(0)-T4(r) 



(4.2.39) 
(4.2.40) 

(4.2.41) 
(4.2.42) 

(4.2.43) 
(4.2.44) 



These eigenvalue equations are much easier to solve than the original Schrodinger 
equations ( 4.2.26| ) and (4.2.32), because the complicated "crystal potential" in 
H m f are absorbed in the effective masses m e , and the bottom and top energies 
of the bands, e c (0) and e v (0), respectively. Once the envelope functions are thus 
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obtained, the wavefunctions can be obtained from (4.2.28) or (4.2.33). Because 
of this simplicity, the Schrodinger eqs. for the envelope functions are widely used 
in the studies of nanostructures. However, care should be taken in discussions of 
optical properties, as we will explain in the subsequent subsections. 

We next consider a state in which an electron-hole pair is excited. In this case 
we must take account of the strong electron-hole interaction. If the band mixing 
caused by the interaction is negligible a pair state (exciton state), which is labeled 
by a set of quantum numbers rj, can be written as 



\xrj) = d 3 r e d 3 r h ip xr) (r e , r h )tpl(r e )^l(r h )\g) , 



<p xr ,(r e ,r h ) = 2^ F xri (k e ,k h )4> ckc (r e )(j)* vkh (r h ) 



(4.2.45) 



(4.2.46) 



Y / J2 F ^( B -o^W( r e - Rj)w* v (r h - Kt), (4.2.47) 



where 

F xv (r e , v h ) = i^'^V k « r «- ik ^F^(k e ,k ft ). (4.2.48) 

The probability amplitude of finding an electron at the position r e and a hole at 
Vh is evaluated as 

(g\tph{r h )tp e (r e )\(p xri ) = (p xv (r e ,r h ). (4.2.49) 

Therefore, <p xr .(r e , r^) can be interpreted as the wavefunction of the pair. It is 
normalized as 



1 = 



(xr)\xrj) = / d 3 r e / d 3 r h \<p xr) (r e , r h )\ 2 = V' V'|F X7) (k e , k h )\ 2 
Jv Jv k e k h 

« 2 VV|F m) (R 3 ,R £ )| 2 = / d 3 r e [d 3 r h \F xv (r e , r h )\ 2 . (4.2.50) 
o Jv Jv 



Note that tp xri (r e ,rh) is not required to be antisymmetric under interchange of 
r e and r^. That is, we can treat an electron and a hole as if they were different 
particles (as far as we can neglect the band mixing). 

The forms ( 4.2.46 )-( 4.2.47 ) of the wavefunction can describe both the Wannier 
exciton (for which the effective Bohr radius a* B 3> lattice constant a) and the 
Frenkel exciton (for which a* B ~ a). We here consider the case of the Wannier 
exciton, and assume, for simplicity, the parabolic dispersions ( 4.2.41 )-( 4.2.42 ). 
Then, in a perfect crystal, the total crystal momentum K of the pair becomes 
a good quantum number, and we write r\ = (K,f), where v denotes the set of 
the other quantum numbers. Considering that K should be shared between the 
electron and hole in proportion to their masses, we can take F xri in the following 
form: 



F xv (k e ,k h ) = X)U(k)« k ., k+ ™ fK «k h) k-^K = Fk»(KM), (4.2.51) 
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where U v is a function to be determined, and M = m e + m/j. Inserting this form 
into Eq. (|4. 2.481), we obtain 



1 



F xr] (r e ,r h ) = ^=exp[iK • R] U v (r e - v h ) = F K „(r e ,r h ) 



where we have introduced the center-of-mass coordinate, 



and 



R 



U v {v) 



m P r„ 



m h r h 



AI 



1 w 



r ^(k). 



(4.2.52) 



(4.2.53) 



(4.2.54) 



We can interpret Eq. (4.2.52) as the product of the center-of-mass envelope func- 
tion e' K ' R /VV and the relative- motion envelope function XJ V . From ( 4.2.50D we 
see that each of them are normalized to 1 : 



exp[iK ■ R] 



/d 3 r|C/,(r)| 2 =^V,(k)| 2 . (4.2.55) 
Jv 



The functional form of t/„ depends crucially on the electron-hole interaction. 
The form of the e-h interaction was reviewed, for example, in 13]. It was shown 
there that the form is nontrivial for r <a* B because of various exchange potentials 
and the finite size of the electron and hole. Since we are considering Wannier exci- 
tons, the long-range parts of the interactions are more important than the short- 
range parts. The dominant e-h interaction would then be the screened Coulomb 
attraction — e 2 /er, where e is the dielectric constant. When the e-h interaction is 
approximated by it, the "Shrodinger equation" for U v takes the familiar form; 



2/y 



er 



U v — E U U U , 



(4.2.56) 



where \x is the reduced mass, = l/m e + l/m/j. By solving the simple eigenvalue 
equation (4.2.56), we can calculate the eigen-energy Evl v and the wavefunction 
<f¥Lv of the state |Kj/). That is, 



Ekv = Eg + E v 



h 2 K 2 

2M ' 



(4.2.57) 



We see that the exciton energy consists of the band gap energy Eq, the energy of 
the e-h relative motion E Vl and the energy of the center-of-mass motion h 2 K 2 /2M . 
The binding energy of the e-h pair is —E u . We can obtain the wavefunction 
( Pxri(' r ei r h) = fKi/(r e ,r/,) by inserting the eigenfunction [/„ into Eq. (4.2.52), and 
the resulting expression into ( 4.2.47] ): 



(PKu(r e ,r h ) = v ^2 XI ~7f exp 



/K 



m e R 7 - + TO/iR^ 



AI 



U v {R.j - 



x w c (r e - Rj)wl(r h - TL e ). 



(4.2.58) 
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Or, we may insert U v into ( 4.2.46| ) to get 



V>Ku(r e) r h ) = — !=exp[jK-R] S Kiy (r e ,r h ), (4.2.59) 
v V 



1 v^' 

S-K V {v e ,Y h ) = —f=2-^ exp[ik- (r e -r h )] U v (k) w ck+ ^ K (r e )u* k _m^ K (r h ). 
v v k 

(4.2.60) 

Note that we can replace u c (r e ) and u v (rh) with u c (p e ) and u v (p h ), respectively, 
where p e and denote the position vectors which are reduced in a unit cell (i.e., 
they are r e and measured from the center of the unit cells in which r e and r/, 
are located, respectively). 

Actually, except for high-quality crystals kept at a low-enough temperature, 
the crystal momentum K of the center-of-mass motion is not conserved because 
the exciton is scattered by phonons, impurities and/or defects, and the exciton 
wavefunctions are modified. For example, when the total potential of impurities is 
slowly-varying (i.e., characteristic length of the variation 3> a* B ) the wavefunction 
of an exciton state \Nv) which is trapped by the impurities may be obtained as 
superposition of ^Ki/(re,rfc) over various K's, 

ip Nv (v e ,v h ) = ^2 G N v(K)ip K „(r e ,r h ) (4.2.61) 

K 

\ - \ - , frrieRj+mhRA 

j £ ' 
x w c (r e - Rj)w* v (r h - R e ), (4.2.62) 

where Gnu is a superposition coefficient, with iV being a quantum number, and 

G Nv {R) = -J=^ '^,(K)e iKR (4.2.63) 
\V K 

can be regarded as the envelope function of the trapped center-of-mass motion. It 
is normalized as 

/ G* Nv (R)G N > v (R)d 3 R = Y^G* Nv {K)G N . v {K) = S N , N >. (4.2.64) 
Jv K 

The trapping may be viewed as the Anderson localization of the exciton, where 
the size £g of the spatial extension of G^r^R) is the localization length £i oc . We 
may also use the form fl4.2.6l| ) or ( [4.2.62] ) as an approximate wavefunction when 
the quantum coherence of the center-of-mass motion is broken by scatterings by 
phonons or other excitons. In this case £q is the phase breaking length Both 
the localization and the phase breaking may be important in real systems, and Iq 
would be given roughly by 

£ G ~min[^ oc , (4.2.65) 

It will be shown later that the oscillator strength of the exciton is proportional to 
£q when £q is less than the wavelength of the optical field. 
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In a similar manner, states with two or more e-h pairs can be written in terms 
of two- or more-particle wavefunctions as 

/ d 3 r\ [ d\\ ■■■( dPri [ d\l ■ ■ ■ rtrl,4, - ,'1^,-) 
Jv Jv Jv Jv 

X V-t (r l)^t (r 2 ) ...^t (r l ) ^t (r 2 ) ...| 5) _ (42 g 6) 

General states are linear combinations and/or classical mixture of such states with 
zero, one, and more e-h pairs. 

4.2.3 Optical Matrix Elements 

The envelope functions should not be taken as the true wavefunctions because 
optical phenomena arc frequently accompanied with transitions between the cell- 
periodic parts u&k's of the Bloch functions. That is, optical matrix elements depend 
not only on the envelope functions but also on the functional forms of Ubk's (or, 
equivalently, of the Wannier functions). The latter dependence becomes crucial 
when nonlinear optical processes or intraband transitions are discussed. Unfor- 
tunately, this point was disregarded in much literature, and unphysical results 
were sometimes reported. We here present typical matrix elements for the states 
described in the previous subsection. 

The interaction with an external optical field A(r, t) in the Coulomb gauge 
takes the following form: 



Ui = fd 3 rft(r) 
Jv 



2 



■- — A • V + - — 

imc 2mc z 



^(r), (4.2.67) 



= n { p + nf\ (4.2.68) 

where TL^p is linear in A, whereas Tif^ is quadratic. To evaluate various matrix 
elements of f6p and Hf \ it is useful to derive formulas for matrix elements of a 
general operator of the form 

Q= [d 3 r [ d 3 r' Q(r,r')^(r)ij(r'). (4.2.69) 
Jv Jv 

In some cases, its expectation value in the ground state \g) does not vanish. Such 
a non-zero expectation value in the ground state can be absorbed by redefinition 
of the operator. This means a change of the zero of energy when Q is a part of the 
Hamiltonian. The redefinition is most conveniently accomplished by the "normal 
order" with respect to i/> e and iph- The normally-ordered Q is 

: Q : = [d 3 r f d 3 r' Q(r, r') : ^(r)^!-') : 
Jv Jv 

= [d 3 r [ d 3 r' Q(r,r') Ul{v)^ & {v') + $( r )$ (r ') 
Jv Jv L 

+^(r)4(r') - ^(r')^(r)l , (4.2.70) 
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which differs from Q only by a c-number term. Using Eq. (4.2.14) and the anti- 
commutation relations, we can easily show that 



(g\-Q- Is) = o, 



(xrj\ :Q:\g)= dn dr 2 tp* x (ri,r 2 )Q(ri,r 2 ), 
Jv Jv 



(4.2.71) 
(4.2.72) 



(xri\ : Q : \xrf) = / d 3 n / d 3 r 2 / d 3 r 3 [ip* x (r 1 ,r 3 )Q(ri,r 2 )^/(r 2! r 3 ) 
Jy Jv Jv 

-<P*x V ( r 3, ri)Q(r 2 , ri)^(r 3 , r 2 )]. (4.2.73) 

These formulas enable us to evaluate various matrix elements from the two-body 
wavefunction <p xv (r e , r/J. 

For an optical field whose photon energy is of the order of Eg or less, the 
wavelength A of A is hundreds nano meters or longer. When the spatial extension 
of ifxrii^e, f/i) is smaller than A we can approximate matrix elements of : TLi : as 

(x V \ : H { P : \g) ~ -— A(r x „) / d 3 r (x V \ :ft(r)e-V $(r) : \g), (4.2.74) 

2 

(xr/l : ttf : \g) ~ -^A 2 (r„ ; ) / d 3 r : ^(r$(r) : | 3 >, (4.2.75) 
2m^ 7 y 

and similarly for other matrix elements. Here, r xr/ denotes the center position 
around which tp xri is localized. We may use, for example, ifNv for tp XT ], and evaluate 
the above integrals using (4.2.72) by putting Q(r,r') = e-V<$(r — r') and Q(r,r') = 
5(r - r'), respectively. From Eqs. Q4.2.59D , fl4.2.60| ) and ( |4.2.6l| ), we find 



f d 3 r (Nu\ : ft(r) e • V </>(r) : \g) = / ' d 3 r [ e ■ V ^(r,r')] r ,^ r (4.2.76) 
Jv Jv 

= G Nu (0)J2'u:(k)(ck\ e ■ V |«k), (4.2.77) 



d 3 r (Nv\ : ^(r)^(r) : |g) = / d 3 r ^„(r,r) = 0, 
where we have utilized Eqs. fl4.2.3| ), ( [4.2. 5| ) and ( |4.2.38| ), and 

(ckl e • V luk) = — / d 3 p u* k (r) e ■ V u„ k (r) 

"0 Jvn 



(4.2.78) 



(4.2.79) 



represents contribution from each cell. If this integral is independent of k in the 
range of k for which f/„(k) is not small, then we denote the integral by (c| e • V \v) 
and Eq. (4.2.77) is reduced to 



d 3 r (Nu\ : ft(r) e . V ^(r) : \g) = W G^(0)[/*(0)(c| e • V \v) (4.2.80) 
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At first sight this result might look unphysical because of the anomalous factor 
VV '. However, this factor is canceled by the factor 1/vF in Gat„(0). This can be 
seen as follows. Recall that we have assumed that the exciton state is well localized, 
Ig < A. For ipnl, to be well localized Gnv(K) must be a smooth function of K, 



and we can replace the summation in Eq. (4.2.64) with the integral as 

1 = J2'\Gn„(K)\ 2 = ^3 J d 3 K \G N „(K)\ 2 (4.2.81) 

Considering also that the spatial extension of Gnu is of the order of £q, we find 
that the maximum value of |G/vv(0)| should behave as 

max 1^^(0)1 ~ Jy- (4.2.82) 

where we have put, apart from an unimportant numerical factor of order unity, 

V G ~ £% (4.2.83) 



Therefore, the matrix element (1.2.8C) is independent of V, as it should be. The 
matrix element is actually proportional to \JVg, which is sometimes called the 
"coherence volume" [|| || because it gives a measure of the volume of the space 
region over which the exciton wavefunction is coherent. An important consequence 
of the Vg dependence of the matrix element is that the radiative recombination 
rate T 1 ^ of the exciton is proportional to V G in the region where £g < A 



Nv 

According to Eq. ( 4.2.6E| ), this phenomenon may be observed in quantum 

^rad 



wells and in quantum dots, as the dependence of T 1 ^ on temperature (on which 
depends) and on the dot size (~ £ioc), respectively. The size dependence of 
of quantum dots was clearly observed by Nakamura et al. || and Itoh et al. 
0, whereas the temperature dependence of of quantum wells reported by 
Feldmann et al. [Q could be interpreted as due to the temperature dependence of 
the exciton distribution because the experiment was performed at relatively high 
temperatures 



pendence of u c and u* in Eq. ( 4.2.59 ) as 



Note that the result ( 4.2.8C) can also be obtained when we neglect the k de- 



S Ku (r e ,r h ) ~ U v (r e - r h )u c0 (r e )u* v0 (r h ) = 5^ 0) . (often wrong) (4.2.84) 

In this approximation, the wavefunction ifNv, for example, is decoupled into the 
envelope function of the center-of-mass motion, that of relative motion, and cell- 
periodic parts as 

(pNv(r e ,r h ) ~ G , Ar ly (R)f/ I ,(r e - r h )u c0 (Y e )u* v0 (v h ) = <p$ v (r e ,r h ). (often wrong) 

(4.2.85) 

This approximation has been used in much literature because of its simplicity. 
This approximation is sometimes good for simple cases such as the band-to-band 
matrix element between direct-allowed bands which are non-degenerate. In studies 
of nonlinear optical phenomena, however, we must also evaluate other matrix 
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elements. For such general matrix elements, the approximate forms (1.2.84) and 
(1.2.85) often give wrong and unphysical results, as we now explain. 

As a typical example let us investigate the transition matrix elements between 
different exciton states. When the spatial extension of both the initial and final 



states are smaller than A, we find from formula (4.2.73) 



(xr,\ : H [ P : \xr,') ~ -^±A(v xri ) / d 3 r e f 
T-mc J v J v 



d 3 r h ^>x V ( r e,r h )e-(V e +\7 h )ip xri ,(r e ,r h ). 

(4.2.86) 

Here, you can arbitrarily replace A(r XJ1 ) with A(r xv '), because A(r XJ1 ) ~ A(r xv ') 
from our assumption A < Iq. The point is that derivatives of both the envelope 
functions and cell-periodic parts contribute to this matrix element. For the latter, 
we can use the group-velocity formula; 



d d ^ fc *k(r) — V0 6k (r) = - 



1 de b (k) 



h <9k 



(b = c,v) 



(4.2.87) 



where m denotes the free-electron mass. For the parabolic dispersion (4.2.41) 
fl4.2.42| ), this yields 



(ck|V|ck) 
(vk|V|uk) 



--l)k, 



m 



m h 



l k. 



(4.2.88) 
(4.2.89) 



It is seen that the derivatives of the cell-periodic parts can never be neglected 
unless m e = m and uih = —m. For example, when m e ~ m/j ~ m the derivatives 
are (ck|V|ck) ~ and (wk|V|wk) ~ — 2ik. Physically, these derivatives represent 
modification of velocities (whose operator is hV jim in the absence of electromag- 
netic fields) due to Bragg reflections by the crystal potential. The modification is 
larger for holes than for electrons because usually both m e and are positive. 
Use of the simplified forms (4.2.84) or (4.2.85) leads to the total neglect of these 
derivatives (because they vanish at k = 0), and often results in unphysical results. 
For example, when (p xv — tpNv we find 



(Nu\ : H { P : \N'v') 



imc 



v y f d 3 i?G^(R) e- V GW(R) 
Jv 



-- I d 3 R G% l/ (R)Gjsr' V '(R) fd : 
M Jv Jv 



r tC(r) e- V U v ,{r) 



.(4.2.90) 



Note that the first and second terms of the right hand side (rhs) become finite for 
different values of N, v, N' and v' . Note also that the typical value of the first term 
is much smaller than that of the second term because a* B <C Ig (by our assumption 
on ipNv) and fi < m (e.g., m/fj, ~ 10 for GaAs). On the other hand, if we took 
^xr) = ¥>a/v we could obtain the first term only; 

(Nv<®\ : H ( P : |JVV(°)) ~ A{r xn )5 u y [ d 3 R G* Nu (R) e ■ V G N ,„(K). 

imc J v 

(4.2.91) 
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It is clear that the use of the simplified form would lead to wrong results 
for the matrix element, both in magnitude and in the selection rules. Therefore, 
we cannot use the simplified forms ( 4.2.84 ) and ( 4.2.85 ) for physical processes in 
which exciton scatterings by optical fields are involved. Such processes include, for 
example, second-order nonlinear responses and third-order processes in which the 
exciton scatterings by optical fields are involved as intermediate (virtual) states. 

It seems worthwhile to point out another difficulty arising from the use of the 
simplified forms. Equations ( 4.2.88| ) and (4.2.89) suggest that u&k's play important 
roles in the gauge invariance of the theory. However, the use of the simplified forms 
is equivalent to the fixing of u^k's on itbo's which are defined in some particular 
gauge. As a result, the gauge invariance of such a theory is violated, ^and physical 
predictions of the theory depend on the gauge [jlO| pj| . No definite conclusions 
can be drawn from such a theory, of course. 

Both of the above difficulties become severer if one proceeds from linear opti- 
cal processes to nonlinear processes. Unfortunately, however, the simplified forms 



(4.2.84) and (4.2.85) have been used in much literature which discusses nonlinear 
optical processes, without any consideration on these difficulties. The most appro- 
priate way of calculation is to use the correct wavefunctions ( [4.2.6C ) and (4.2.61) 
together with the original interaction Hamiltoni an (4~2.68| ). There are alternative 
ways in which one may use the simplified form ( 4.2.85 ). One way is to use of the 
"length formula" for the interaction M [fTo|: 



njlength 

it j — C 



d 3 r $+(r) r • E ^(r). 



(4.2.92) 



This expression is rather anomalous because the value of r extends to infinity as 
V — > oo. Therefore, the normally- order form : -f[ l ^ n 9 th • should be used in ac- 
tual calculations. I t is no t difficult to show that this interaction together with 
the simplified form ( 4.2.85| ) give the correct results. The other way is to perform 
calculation according to the following prescription: In calculations of the envelope 
functions and matrix elements of intraband transitions (exciton-exciton scatter- 
ings), use the following Luttinger-Kohn Hamiltonian as an effective "intraband 
Hamiltonian" for the envelope function G^yUy [fLOl; 



rrenv 
intra 



2m P 



fir, 
-V e 

I 



e -A 



2mh 



-V h --A - 
i c ) e|r R 



(4.2.93) 



In calculations of matrix elements of interband transitions (exciton creation and 
annihilation), on the other hand, one must use the bare interaction Tij. Physically, 
this prescription indicates that slowly-varying phenomena are determined by the 
effective masses whereas the rapidly-varying phenomena are determined by the 
free-electron mass [EJ . 

2 There is another possible origin of the violation of the gauge invariance of the theory. It 
is some nonlocal interaction, such as the exchange interaction, which may appear as a result of 
a mean-field approximation [El. In our case, however, such an interaction has been neglected 
because we are considering Wannier excitons. 
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4.2.4 Quantum States in Nanostructures 

The above formulation must be modified to treat nanostructures, such as quantum 
wells (QWs) and quantum wires (QWRs), because the band structures are mod- 
ified by the quantum confinement, and also because nanostructures are usually 
made of semiconductors which have degenerate valence bands (whereas we have 
assumed for simplicity that the bands are non-degenerate). One way to take ac- 
count of these effects is to combine the k • p theory with the exciton theory, which 
would enable us to determine (though perturbatively) both the band structure 
and the exciton wavefunction simultaneously. 

We here describe another simpler way — a rigid-band approximation. This 
approximation may be valid (i) in the case where the quantum confinement is very 
strong (i.e., the confinement energy ^> the exciton binding energy) and also (ii) in 
the case where the quantum confinement is very weak (i.e., the confinement energy 
■C the exciton binding energy) and both the c and v bands are non-degenerate (i.e., 
these bands are energetically separated around k = from other bands, and the 
energy separation is larger than the exciton binding energy). We can take the 
band structures as those of the nanostructure and the bulk crystal for cases (i) 
and (ii), respectively, because the band structure would be rigid against the e-h 
interaction. The electron and hole operators can be defined according to the rigid 
band structure. We can then evaluate exciton states taking the e-h interaction into 
account. Using these exciton states, we can evaluate various optical properties, 
such as nonlinear susceptibilities, of the nanostructure. 

In case (ii), in particular, the rigid bands are the bands of the bulk crystal, and 
the only important effect of the quantum confinement would be the confinement of 



the center-of-mass motion. Hence we may employ ipN„ of (4.2.61) as the exciton 
wavefunction, and the formulations of the previous subsections can be applied 
directly. A typical example is the Is exciton in a quantum dot whose size is 
larger then a* B and which is made of non-degenerate semiconductors. For higher 
energy states (2p, 3d, and so on), the weak-confinement assumption becomes worse 
because the spatial extension of the e-h relative motion becomes larger for higher 
energies. 

In case (i), on the other hand, the band structure would be modified by the 
quantum confinement, but the modified band structure would not be affected 
strongly by the e-h interaction. In GaAs/AlGaAs nanostructures, for example, 
the two-fold degeneracy at the top of the valence bands is lifted by the quantum 
confinement, and two series of subbands are formed, which are called the heavy- 
hole subbands and light-hole subbands. The subbands are no longer degenerate 
at the top of the valence bands. By the assumption that the intersubband energy 
separation is larger than the exciton binding energy, the e-h interaction would not 
alter the subband structure strongly. In this case, we may take the single-body 
eigenstates of the mean-field Hamiltonian of the nanostructure as basis functions. 
For concreteness, we assume for a nanostructure a quantum well structure of well 
width L and area S, and a <C L <C a* B (a is the lattice constant). The cases of 
a quantum wire and quantum dot can be formulated in similar manners. Then, 
a single-body eigenstate (of low energies) may be approximated by the following 
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form; 



'""kn ( r ) = \l-F^2 ex Pi ik \\ 



R^|]F 6 „(^H(r-R^R^), 



(4.2.94) 



where n denotes a quantum number, kn and R^ii denote two-dimensional vectors 



which represent the wavenumber and the cell position, respectively, in the QW 
plane (xy plane), and Zg is the z coordinate of the cell position, i.e., R^ = (R^ , Zg). 
Note that b here labels the bands which are modified (the degeneracy is lifted) by 
the quantum confinement, and Wb{r — R^;R^) denotes the Wannier function of 
the modified band in the well (when R^ G well region) or barrier (when R^ € 
barrier region). The Wannier functions and the envelope function Fb n may be 
obtained, e.g., by using the k- p perturbations. In particular, for bands which are 
non-degenerate in the bulk crystal, we may approximate Wb by that of the bulk 
crystal, and Fb n may be obtained by solving the effective-mass equation; 



£b ( k||, 



d 

' l a~' z 
az 



Fbnk u (z) = ebnk,,Fbnk u (z) , 



(4.2.95) 



where £b(k||, k z ; z) is the dispersion of the bulk band of the semiconductor which 
composes the well (when z € well region) or barrier (when z G barrier region). For 
a parabolic band, in particular, Eq. (4.2.95) becomes 



if 



2?7j fc ||(z) 2uib z (z) dz 2 



e b (0;z) 



Fbnk« (z) = £fcnk,i Fbrik» (z)- 



(4.2.96) 



Here we have considered the fact that in general the effective mass mj, becomes 
anisotropic by the band mixing. The k|| — > limit of this equation is the familiar 
form; 



e b (0;z) 



F bn (z) = e bn Fbn(z), 



(4.2.97) 



2m bz {z) dz 2 

where F& n = Fb n o and £&„ = Sb n o- Equations ( |4.2.96 ) and ( 4.2.97 ) are basically 
particle-in-a-box problems. Note however that the boundary conditions on Fbnk 
or Fbn at heterojunctions depend on the functional forms of Wb{r — R^; R^). For 
example, when Wb takes different values across a hetero junction (because the semi- 
conductors are different on both sides) then Fb n k must be discontinuous in order 
for the wavefunction to be continuous across the heterojunction Jl2[ . 

For a typical GaAs/AlGaAs QW, the conduction band is non-degenerate and 
we may use Eq. (4.2.95), whereas the valence bands are degenerate and we should 
use another method such as the k-p perturbation. However, Eqs. (4.2.95)-(4.2.97) 
have been widely used even for the valence bands. This may be justified if we 
replace Wb with that of the QW (i.e., Wb which is evaluated with the k • p per- 
turbation) and simultaneously adjust parameters in £b(k|| ,k z ;z) of Eq. (4.2.95) in 



such a way that the wavefunction agrees with the result of the k • p perturbation. 

We can also construct the wavefunction of an exciton, using the Wannier func- 
tions of the QW, as 



ip xv (r e ,r h ) 



v ^2^2 F x „(Rj ,R e )w c (r e - RjiR^w^Th 

3 t 



R*;R*). (4.2.98) 
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The envelope function F xri (r e , r^) can be evaluated from, say, an effective-mass 
equation. Using (/J^'s thus obtained, we can evaluate various optical matrix ele- 
ments in a similar manner to the previous subsections. The most important point 
in this case is that the Wannier functions (and therefore the cell periodic functions) 
become anisotropic when the band degeneracy is lifted by the quantum confine- 
ment p^ |. For example, when the valence band degeneracy is lifted, the matrix 
elements between the cell periodic functions become anisotropic, and, for small k, 

(771 \ 
— -l k, (4.2.99) 
m h (k) J 

where k denotes the unit vector which is parallel to k. The interband matrix 
element also becomes anisotropic, which for small k we denote as 

lim(ck|e • V|uk) = (c\e ■ V|uk). (4.2.100) 

These anisotropies of the cell-periodic parts have been observed as anisotropics of 
various optical properties of QWs. It is clear that as is the case of the bulk crystal 



(subsection 4.2.3) optical matrix elements must be evaluated carefully. 

We have calculated in Refs. Q [|l5] Q Q Q linear as well as nonlinear 
optical properties of nanostructures using the above formulation. Up to now some 
of the results have been confirmed experimentally [^0| Jl{j] |2lJ . 

4.2.5 Quantum Optical Phenomena in Nanostructures 

So far the photon field has been treated as a classical field. In recent years, 
however, quantum optical phenomena in nanostructures have been studied by 
many researchers. We here depict some examples from our work. 

The first example is the quantum nondemolition measurement of the photon 
number using an electron interferometer which is composed of double quantum 
wires [^2| [ p3| . It was shown that one can measure the photon number of a quantum 
state of photons using the electron interferometer as a photodetector in such a way 
that the post-measurement state of the photons carries the same information on 
the photon number as the pre-measurement sate. When the pre-measurement 
state is an eigenstate of the photon number, in particular, the wavefunction is 
unchanged by the measurement, whereas the observer obtains information on the 
photon number. In contrast, when the observer uses conventional photodetector 
the information on the photon number is lost by the measurement: when he uses 
a photodiode, for example, the post-measurement photon state is the zero-photon 
state, which has no information on the pre-measurement state. In Refs.[^2| p3)] 
various quantities, such as the density operator of the post measurement state and 
the measurement error, have been calculated as a function of the pre-measurement 
state of the quantized photon field and the structural parameters of the quantum- 
wire interferometer. 

The second example is the photon creation from a false vacuum of semiconduc- 
tors p4[ . In quantum field theory, particles are created when a parameter(s) in the 
Lagrangian (or, equivalently, in the field equations) have time dependence even if 
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the initial state is the vacuum state (which is defined before the time variation 
begins) p5| . Since this is a general property, we can expect such phenomena also 
in condensed matter. For example, photons may be created if a material parame- 
ter (s) varies as a function of time (2(| | p7[ . We have pointed out that conventional 
discussions on the photon creation in condensed matter, which are based on phe- 
nomenological quantization of the macroscopic photon field, break down near the 
exciton resonance [^4| . To treat the exciton resonance correctly, we have proposed 
the use of a two-field model, in which both the photon field and the exciton field 
is taken into account as microscopic degrees of freedom . The two-field model 
has allowed us to evaluate the creation spectra over a wide energy range. It is also 
pointed out there that even a slight singularity in the functional form of the time 
variation of the parameters causes large deviation from non-singular (physical) 
results, whereas singular forms were assumed in much literature. 

The third example is the radiative lifetime of an excited atom in an inhomoge- 
neous and absorptive cavity ||^] . We have demonstrated strong dependence of the 
lifetime on the material parameters such as the absorption coefficient and the cav- 
ity structure. The forth example is the photon statistics of light emitting diodes 
(LEDs) at a low injection level [^9|. We have derived a formula which gives the 
photon statistics as a function of the pump statistics, device parameters, and the 
detection efficiency of each photon mode. 

These are only a very limited examples of the quantum optical phenomena in 
nanostructures. More examples will be presented in the subsequent sections and 
in the references cited therein. 
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